Methode van Roberts - hidden line/hidden surface 
algoritme 

Computergrafiek - Computer graphics 

Inleiding 

De methode van Roberts is het eerste hidden line algoritme dat gemaakt werd. Het algoritme 
werd in 1963 door Roberts “uitgevonden”. 


Doel en veronderstellingen 

De bedoeling is om in een 3D ruimte te bepalen wat zichtbaar is en wat niét zichtbaar is. We 
veronderstellen N voorwerpen in de ruimte die elkaar niet snijden en die eveneens allemaal 
convexe hoeken hebben, (met andere woorden: geen enkele hoek mag > 180° zijn. Indien dit 
wel het geval is, kan je met een ander algoritme de ruimte eerst verdelen in kleinere 
voorwerpen, zodat de niet-convexe voorwerpen (ook wel concave voorwerpen) worden 
gesplitst en zo allemaal wél convexe voorwerpen worden). 

Het Robert algoritme gaat er eveneens van uit dat al onze voorwerpen niét doorzichtig zijn. 
Tot slot gaan we ervan uit dat de voorwerpen bestaan uit veelvlakken, bollen of cirkels laten 
we dus ook buiten beschouwing. (Ter info: je kan een cirkel of bol eventueel gaan voorstellen 
door zéér veel kleine driehoekjes/vierkantjes, die zo visueel hetzelfde voorwerp lijken voor te 
stellen) 

We willen dus onze 3D ruimte gaan voorstellen op een 2D ruimte, met andere woorden een 
3D ruimte op een blad of een computerscherm gaan voorstellen. 

We gaan dus voor alle ribben van de voorwerpen gaan nazien of deze nu wél of niét zichtbaar 
zijn. (Een ribbe is de plaats waar 2 vlakken elkaar snijden.) 
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Doel is dus te weten dat op bovenstaande figuur het rood aangeduide wél zichtbaar is, maar 
een stukje niét zichtbaar is. 

Onzichtbaar? 

Een ribbe is onzichtbaar in 2 gevallen: 

• Ze wordt bedekt door het voorwerp (veelvlak) waar ze zelf deel uit maakt, (denk aan 
een doos die schuin voor je staat, de achterkant kan je dan niet zien, aangezien de 
voorkant je het zicht belemmert. 

• Ze wordt gedeeltelijk of geheel bedekt door een ander voorwerp in de ruimte. 

Het is mogelijk dat een lijnstuk in vele verschillende stu kk en zo visueel ‘geknipt’ 
wordt (ze zeggen dat het niet noodzakelijk enkelvoudig samenhangend is). 

Denk maar aan een borstel die voor een doos staat. Tientallen/honderden dunne 
balkjes (de haren van de borstel) delen de ribbe van de doos in tientallen/honderden 
stukjes. 

We willen dus bepalen of een rechte (ribbe) nu wél of niet zichtbaar is, en dat dankzij de 
methode van Roberts. 


Enkele definities 

Alles laten we gebeuren in de zogenaamde co-ruim te (“omega ruimte” uitgesproken). 

[hiervoor heb je al wat grafische voorkennis voor nodig. De hele ruimte is bepaald in een 
gewoon XYZ-assenstelsel. We nemen vervolgens in die ruimte een oogpunt (punt waar we 
gaan kijken in de ruimte), E genoemd, en VOOR ons oog (punt E) zetten we loodrecht een 
scherm. Dat scherm, een vierhoek, is de «-ruimte. Het is dus het scherm waarop de hele 3D 
ruimte zal worden afgebeeld, die gezien wordt vanuit ons standpunt (het oogpunt E)] 

Een punt P in onze ruimte heeft de coördinaten: 

P = (XW, YW, ZW, W) 

Dit zijn zogenaamde homogene coördinaten. 

[Een punt (X,Y,Z) is voor te stellen in homogene coördinaten door gewoon (X, Y, Z, 1) van te 
maken.] 

Een vlak n (pi) in onze ruimte wordt bepaald door 4 coëfficiënten en geschreven in een 
matrix: 


7T = [ a ] 

I b I 
I c I 
[d] 

Verder hebben we nog een hele leuke eigenschap. Een willekeurig punt P in de ruimte ligt 
OP ons vlak n indien het scalair product van P en n gelijk is aan 0. 

Dus: 

P . n = 0 
Of anders geschreven: 
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(XW, YW, ZW, W) . [ a ] = O 

I b I 
I c I 

[d] 

En deze matrixvermenigvuldiging, wordt als volgt uitgerekend en het volgende is er dus 
equivalent met: 

aXW + bYW + cXW + dW = 0 

Al het voorgaande delen door W (linker en rechterlid delen door eenzelfde van 0 verschillend 
getal) geeft: 

aX + bY + cZ + d = 0 

Dus indien aan voorgaande vergelijking voldaan, dan ligt het punt OP n, anders niet. 

Indien ons punt niet op n ligt, dan ligt het er dus boven of onder. Dit maakt dus dat 
voorgaande formule, dan zal zijn: 

P . 7T > 0 of P . 7T < 0 

Verder zegt men dat het vlak n “gericht” is. Het vlak (bekijk het als een gigantisch oneindig 
groot blad papier dat de ruimte in 2 deelt), heeft een boven en een onderkant. Men noteert dit 
als een positieve en een negatieve kant. 

Men kan het vlak omdraaien (dus positief en negatief omwisselen) door de coëfficiënten van 
het vlak te vermenigvuldigen met een bepaalde k < 0. 


Volumematrix 

We gaan verder met nog enkele definities te maken voordat we aan het algoritme zelf kunnen 
starten. 

We definiëren een veelvlak Bi (dus een voorwerp in de ruimte) bestaande uit n, vlakken. 

Dit maakt dat: 

Bi — (7Til TCi 1 TCi3 ••• TCn) 

We hebben dus N voorwerpen in onze ruimte, i gaat dus van 1 tot N (ingewikkelder 
genoteerd: Bi (Vi = 1, ..., N) ) 

We gaan ons veelvlak (voorwerp) zodanig gaan definiëren dat alle vla kk en waaruit het 
bestaat, hun positieve kant naar BINNEN gericht zijn. Indien dit niet het geval is, wordt het 
vlak dus vermenigvuldigd met een k < 0 zodat het wel goed gericht is. 

Dus IN ons voorwerp is het één en al positief, en BUITEN ons voorwerp is het langs alle 
kanten negatief. Dit is belangrijk, want het hele Roberts algoritme steunt hierop. 

Verder hebben we ook nog onze zogenaamde “kijkbalk”. Dit is de balk die we zien. We zien 
namelijk op ons vlak en dat heeft een linkerkant, rechterkant, bovenkant, onderkant en een 
voorkant. Onze kijkbalk is half oneindig; vertrekkend vanuit E tot in het oneindige, [de 
benaming “halfoneindig” lijkt wel vreemd, maar dat is dus dat het ergens start, en dan in het 
oneindig verder gaat. Helemaal oneindig zou zijn dat het in het oneindige start, en tot in het 
oneindige doorloopt; wat dus duidelijk verschillend is met halfoneindig] 
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Met andere woorden: 5 vlakken. En deze 5 vlakken worden opgegeven in onze co-ruimte: 
Xco = 1 , Xco = -1 

Yco = 1 , Yco = -1 

Zco = 0 


We hebben de volumematrix Bo van de kijkpiramide. Dit is dus de matrixvoorstelling van de 
5 vlakken die onze kijkbalk voorstellen: 



0 o 

-1 1 

o o 

1 1 


o 

o 

-1 

o 






Verder gaan we de oriëntering van onze kijkbalk na of deze echt wel goed is (of dus alle 5 de 
vla kk en wel met hun positieve kant naar binnen liggen) door het controleren dat onze Z-as (de 
richting waar we naar kijken) aan de positieve kant ligt. 


Een bijzonder belangrijke eigenschap voor onze methode van Roberts is de volgende. 

Een punt P is zichtbaar indien het IN de kijkpiramide ligt, en dus onzichtbaar als het buiten de 
kijkpiramide ligt. Dit controleren is doordat het punt aan de positieve zijde moet liggen van 
alle 5 vlakken van de kijkpiramide. 

Denk hier even goed over na. Alle 5 vla kk en zijn met hun positieve kant naar elkaar gericht. 
Als ons punt dus voor alle 5 de vla kk en langs de positieve zijde ligt, ligt het dus IN de 
kijkpiramide, en is het dus zichtbaar. 


We kunnen het voorgaande ook noteren als volgt: 
( P . Bo) > 0 


Dus de coördinaten van P vermenigvuldigen met de matrix van Bo, en indien die uitkomst 
groter is dan 0, ligt het dus aan de positieve kant van alle 5 de vla kk en van onze kijkpiramide. 


Het echte werk 

We gaan nu met het echte werk beginnen voor de methode van Roberts. 

We nemen een rechte in onze ruimte: P 1 P 2 . We nemen op deze rechte een generisch punt P: 
P = ZP1 + (1 -X) P2 
met X in [0, 1] 

Met andere woorden: het punt P ligt TUSSEN Pi en P 2 en de positie waartussen het ligt, 
wordt bepaald door een parameter X die tussen 0 en 1 ligt. 

Dat punt P op onze rechte is nu ZICHTBAAR voor ons indien de coördinaten van P 
vermenigvuldigd met elk van de 5 vla kk en van onze kijkpiramide, telkens > 0 geeft. Dus dat 
P langs de positieve kant ligt van alle 5 de vlakken van onze kijkpiramide. 

In formulevorm geschreven: 
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(XPi + (1-A,)P 2 ) . 7t 0l > O , Vi= 1, ... 5 

Dus: P . 7toi moet groter dan 0 zijn, en dat voor alle 5 de vlakken n van Bo wat onze 
kijkpiramide is. 

We herschrijven dit: 

X ( (Pl-P 2 ) . TlOi ) > - ( P 2 . TlOi ) , Vi = 1, ... 5 


En indien we dit opnieuw herschrijven door het vervangen van enkele uitdrukkingen door 
nieuwe variabelen: 

Stel: 

Pi = (P 1 -P 2 ) • 7t0i 
ti = - ( P 2 . 7t0i ) 

Dan: 

Api > ti , Vi = 1, ... 5 

Afhankelijk van het teken van pi is dit een boven of een ondergrens voor X, 

Hier moet je ook eens even goed over nadenken om dit te begrijpen. 

Indien pi > 0, en we verplaatsen die naar de rechterkant van onze vergelijking, dan zegt die: 
X > (h)/ (pi) , Vi = 1, ... 5 

En geeft die dus een ondergrens, want die zegt dat X groter is dan de rechtse uitdrukking. 

Indien pi < 0, en we willen die naar het rechterlid verschuiven, dan dienen we de > te 
vervangen door een < (dat is gewoon zo volgens de wiskunde), en krijgen we dus volgende 
uitdrukking: 

X < (ti)/ (pO , Vi = 1, ... 5; en (p; < 0) 

En geeft die uitdrukking dus een bovengrens aan, want de uitkomst van de rechterkant is 
groter dan X en is dus X kleinder dan het rechterlid. 

Op deze manier, door het gebruik maken van deze formules, bepalen we zo elk gedeelte van 
een ribbe met een functie, ook wel “parametrisch bepaald” genoemd. 


Ribbe onzichtbaar door voorwerp zelf 

We nemen 2 vla kk en die tegen elkaar liggen en die dus een doorsnede (intersectie) hebben. 
Op die intersectie hebben we een ribbe. Definieer het ene vlak als TXij en het andere als K\k', 
met dus “i” als gemeenschappelijke. 

Beschouwen we nu een ribbe Rijk die dus gevormd wordt door de intersectie van de twee 
vlakken, van het voorwerp Bi. 

[Merk op: niet alle vla kk en vormen 2 aan 2 een ribbe, bijvoorbeeld 2 evenwijdige vlakken 
hebben helemaal geen doorsnede en dus ook geen gemeenschappelijke ribbe]. 

Deze ribbe Rijk is Onzichtbaar achter het voorwerp Bi waarvan het zelf deel uitmaakt indien 
het waarnemingspunt E (dus ons oog) aan de POSITIEVE kant ligt van zowel Ttij als k±. 

Ah ja, natuurlijk! Als het langs beide kanten langs de positieve kant ligt, dan wil dit zeggen 
dat beide vla kk en met hun negatieve kant weg van ons zijn (en de negatieve kant is de 
zichtbare kant), en dus is de ribbe bedekt door het voorwerp zélf. 
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Kijk maar naar een doos. Neem een ribbe die vanuit jou standpunt onzichtbaar is. Beide 
vla kk en waaruit die bestaat, liggen allebei met hun positieve kant naar jou gericht! 

Merk op dat één van beide met hun positieve kant naai'jou niet genoeg is, want dan is het 
andere vlak nog wel met de negatieve kant naar je toe gericht én is het dus zichtbaar (want de 
negatieve kant naar je toegericht is nu net per definitie zichtbaar). 

De volgende figuur verduidelijkt dit. 




Wat je ziet op de figuur zijn dus 2 vlakken, die samen de ribbe Rijk vormen. Ze zijn met hun 
positieve zijde naar het gezichtspunt (E) gericht, en zijn dus onzichtbaar. 

Ribbe onzichtbaar door andere voorwerpen in de ruimte 

We bekijken nu het tweede geval: wanneer een ribbe niet bedekt wordt door het voorwerp zelf 
waarvan het deel uitmaakt, maar dat het wordt gedeeltelijk of geheel onzichtbaar wordt 
gemaakt door andere voorwerpen in onze ruimte. 

We gaan opnieuw een generisch punt P nemen op een lijnstuk P 1 P 2 . 

We gaan ook weer 2 vlakken beschouwen en hier ook weer een ribbe Rp qr op beschouwen; 
helemaal net zoals hiervoor. 

Onze P is dus opnieuw: 


P = kPi + (l-k)P 2 


We willen dus de waarden van k zoeken waarvoor het lijnstuk onzichtbaar is achter een 
voorwerp Bi. Die k geeft dus de positie van P aan op het lijnstuk P 1 P 2 en dus van waar tot 
waar het zichtbaar/onzichtbaar is. (en we gaan dat dan uiteindelijk doen voor alle voorwerpen 
in de hele ruimte, zodat we kunnen beslissen wat er wel/niet zichtbaar is van het lijnstuk). 

We gaan een niéuw punt beschouwen, dat we Q gaan noemen. Dit punt (in homogene 
coördinaten weer), wordt als volg beschouwd: 
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Q = P + a E 
ofwel 

Q= X Pi + (1-1) P 2 + a E met a>0 

en E = (0, 0, -1, 0) (waarnemer op - 00 ) 

Dit vraagt om een beetje extra uitleg. 

Q is dus een functie van P, en eveneens van E (het oogpunt). Punt P hangt dan weer af van X. 
En dat alles waarbij de waarnemer in de homogene coördinaten E = (0, 0, -1, 0) staat, of met 
andere woorden op -00 (min oneindig). 

Voor de volledigheid, de technische definitie zegt: 

“Q ligt op een halfoneindige strook rustend op P1P2 , en beschreven door een rechte die elk 
punt van P1P2 verbindt met E”. 

Een tekening maakt veel duidelijk. 



Wat we zien op de tekening is dus de blauwe pijl die onze kijkrichting aangeeft met het rode 
bolletje, wat ons oogpunt E is. Het blauwe parallellogram is de co-ruim te, met het assenstelsel 
X, t „ Y„, en Z M op aangeduid. De groene stippellijnen geven de kijkpiramide aan. 

We zien dan de rechte P 1 P 2 , met erop een punt P, dat afhankelijk is van 1, vandaar P(l) 
genoteerd. 

Vanuit P(l) vertrekt er een rechte die naar het oogpunt E gaat, en die op een bepaalde plaats 
stopt, en daar ligt het punt Q, dat dus niet enkel afhankelijk van X (want het hangt af van 
P(k)), maar tevens van ot, want a bepaalt de afstand ten opzichte van P(l). 

Q (k, ot) ligt dus altijd op een rechte tussen het oogpunt E en het punt P(k). Het rode/roze is 
het vlak dat de rechte P 1 P 2 en E verbindt met elkaar (en waar Q dus opligt). 

En nu komt weer een idee van Roberts methode: 
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P(A) is ONzichtbaar (het ligt dus achter het voorwerp Bi) indien er een waarde voor a bestaat 
zodat Q een inwendig punt is van Bi. 

Dit vraagt ook weer om wat extra uitleg. 

Wat we eigenlijk doen is het punt P( X) naar ons oogpunt toekomen; en dat doen we dus met 
Q. We laten dat punt dus als het ware steeds dichter en dichter komen, en dat doen we door a 
te laten veranderen. Als we ooit een waarde voor a kunnen vinden waardoor ons punt Q 
BINNENIN het voorwerp Bi ligt, dan is ons punt P(k) onzichtbaar. Natuurlijk, want dan wil 
dit zeggen dat het voorwerp Bi in de weg zit tussen ons oog en het punt P(k). 

Neem weer de doos. En zet nog een ander voorwerp tussen jou en de doos. Neem een punt 
op een van de ribben, en ga in een rechte lijn vanaf dat punt naar je oog. Indien je op een 
bepaald ogenblik tegen dat andere voorwerp komt (en er dus normaal doorheen zou moeten 
gaan), dan wil dit zeggen dat het andere voorwerp dus het zicht blokkeert (in de weg staat) 
van ons gekozen punt en het dus onzichtbaar is! 

We kunnen dit voorstellen met formules. En wel als volgt. 

((XPi + (l-X)P 2 + aE). Ti ij) > 0 
Vj = 1, ... n ; 

Dus het eerste deel van die functie is gewoon Q (want Q = X Pi + (IA) P 2 + a E) en die z’n 
homogene coördinaten vermenigvuldigen we met de coëfficiënten van het vlak Tiij, en dit voor 
alle vla kk en van ons veelvlak (voorwerp) Bi. 

Indien we dus een a hebben gevonden zodat Q langs de positieve kant ligt van élk vlak 
waaruit het voorwerp Bi bestaat, dan ligt het dus IN dat voorwerp en is het dus niet zichtbaar 
voor het oog dat in E staat (want we veronderstelden dat voorwerpen niet doorzichtig zijn!). 

We kunnen die formule nog eens verder herschrijven: 

A. ( (P 1 -P 2 ) . Ttij ) > - ( P 2 • 7tij ) + CC ( E . 7Tij ) 

Vj = 1, ... m 

Dat is dus gewoon enkele parameters van plaats wisselen en wiskundig bewerken, niets 
speciaals aan. 

We gaan nu weer zoals we eerder hebben gedaan, enkele stukken uit deze formule vervangen 
door nieuwe namen, zodat het eenvoudiger lijkt. 

Stel: 

Pij = (P 1 -P 2 ) • 71 ij 
tij = “ ( P2 . TTij ) 

Cij = ( E . TCij ) 

Zodat: 

A pij ^ tij T CC Cij 

Vj = 1, ... n, 
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En ook hier is het weer zo: afhankelijk van het teken van py is dit een boven of een 
ondergrens voor X. 

Indien pij > 0, dan kunnen we het dus naar de andere kant plaatsen en vormt dit: 

X > (hj + ot Cij) / pij 

en is dus X steeds groter dan die uitdrukking rechts en is dit dus een ondergrens. 

Indien pij < 0, dan kunnen we het dus naar de andere kant plaatsen, mits omkering van > naar 
< (is zo volgens de wiskunde), zodat: 

X < (hj + ot Cij) / pij 

en is dus X steeds kleiner dan die uitdrukking rechts en is dit dus een bovengrens. 

Merk op dat de waarde van deze grenzen lineair afhankelijk zijn van a. Maar ten opzichte 
van elkaar wijzigen ze dus niet (door het wijzigen van a wordt een bovengrens niet plots een 
ondergrens of zo). 

Door deze functie krijgen we een hele hoop boven en ondergrenzen. Al deze grenzen bepalen 
het ONzichtbaarheidsgebied voor één welbepaalde a. 


Onzichtbaarheidsinterval 

We definiëren weer eerst enkele zaken. 

We gaan nu de grootste benedengrens nemen, en noemen deze gi,min(a). 

We benoemen ook de kleinste bovengrens, en noemen deze gi,max(a). 

Oi(a) is dan het interval dat afgebakend is door gi, m in(a) en gi,max(a); en dit dus op de “diepte” 
a, want merk op, indien we ot wijzigen (we dus dichter of verder weg gaan van ons oogpunt 
E), kunnen de grenzen en onzichtbaarheid verschillen. 

Een figuur verduidelijkt dit ook weer extra. 



Je ziet dus in het groen de grootste ondergrens, en in het rood de kleinste bovengrens voor X. 
Al de andere zijn andere grenzen die niet interessant zijn. 

Het interval tussen de boven en ondergrens is dan Oi(a), op de figuur in het donkerblauw 
aangeduid. 

De reden dat we dit nemen is goed te begrijpen. We nemen de meest nuttige informatie. 
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Indien we zeggen dat een getal kleiner is dan 10, kleiner is dan 15 en kleiner is dan 7; dan is 
het interessantste om te weten dat het kleiner is dan 7; want die andere 2 gegevens volgen hier 
rechtstreeks uit. 

Indien men zegt dat een getal groter is dan 2, groter is dan 5 en groter is dan 3, dan is het 
interessantste om te weten dat het groter is dan 5, want die andere twee gegevens volgen hier 
ook rechtstreeks uit. 

Zo weetje dan dat het getal groter is dan 5 en kleiner is dan 7; en is ons interval hiertussen. 

We gaan nu Oi definiëren als volgt: 

Oi = Ua>o Oi(a) 

Of met andere woorden: Oi is de unie van alle grenzen van alle a. 

Een tekening kan extra duidelijkheid brengen. 



We zien op deze tekening horizontaal de as voor A, en dat tussen 0 en 1 (dus gaande van links 
naar rechts op ons lijnstuk P 1 P 2 ). 

We zien a en het overeenkomstige Oi(a), wat dus het onzichtbaarheidsinterval is voor die 
welbepaalde a. 

Men heeft nog een tweede uitgetekend, namelijk a’, met ook weer daar zijn eigen Oi(a’) 
onzichtbaarheidsinterval. 

En zo zie je op de tekening steeds links in het groen allemaal puntjes wat de grootste 
benedengrenzen zijn, en rechts in rode puntjes de kleinste bovengrenzen. 

Er zijn dus tientallen boven/ondergrenzen getekenend, enkel halen we er even aena’ uit, 
maar in werkelijkheid doen we dus massaal veel van dergelijke berekeningen. 

Oi is dan de laatste waarde van eender welke Oi(oti) en de hoogste waarde van eender welke 
Oi(oii). En dat zie je bovenaan op de tekening. 

Het is dié waarde die het uiteindelijke onzichtbaarheidsinterval is voor die welbepaalde rechte 
P 1 P 2 . 
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We gaan nu nog gewoon O definiëren tot slot. Deze is namelijk de unie van alle O,: 
O = Ui=l ... N, Oi 


Ah ja, natuurlijk. We hadden eerst Oi(a), dat het onzichtbaarheidsinterval was voor één a. 
Dan maakten we een algemene Oi die het algemene onzichtbaarheidsinterval was voor alle a, 
maar wel ten opzichte van één welbepaald voorwerp B,. 

Nu maken we O, dat het meest algemene onzichtbaarheidsinterval is ten opzichte van alle 
voorwerpen B, met uitzondering van B p , want die moeten we niet meer controleren want onze 
rechte P 1 P 2 , is een ribbe van B p en we hebben eerder al nagezien dat die onze ribbe niet 
onzichtbaar maakt (en dit opnieuw checken zou niet efficiënt zijn). 

Deze procedure moet worden herbegonnen voor élke ribbe R pqr van B p , en voor alle 
veelvlakken B p in onze ruimte. Op deze manier kunnen we de zichtbaarheid van elke ribbe 
kennen ten opzichte van alle andere veelvlakken (voorwerpen) in onze ruimte en hebben we 
het resultaat. 

Dit was het algoritme van Roberts, ofwel de Roberts methode genoemd. 

Bronnen foto’s: 

• PowerPoint slides van Prof. Dr. Ir. Patrick P. Bergmans (Computergrafiek, 
Universiteit Gent, Faculteit Ingenieurswetenschappen). 
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